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Abstract 

The Gauss hypergeometric functions 2-^1 with arbitrary values of parameters are 
reduced to two functions with fixed values of parameters, which differ from the origi- 
nal ones by integers. It is shown that in the case of integer and/or half- integer values 
of parameters there are only three types of algebraically independent Gauss hyperge- 
ometric functions. The e-expansion of functions of one of this type (type F in our 
classification) demands the introduction of new functions related to generalizations of 
elliptic functions. For the five other types of functions the higher-order e-expansion 
up to functions of weight 4 are constructed. The result of the expansion is expressible 
in terms of Nielsen polylogarithms only. The reductions and e— expansion of g-loop 
off-shell propagator diagrams with one massive line and q massless lines and g-loop 
bubble with two-massive lines and q — 1 massless lines are considered. 
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1 Introduction 



The construction of higher order e-expansions of hypergeometric functions has been inten- 
sively discussed in the hterature in the context of the calculation of Feynman diagrams. 
At the present moment, several algorithms for the Laurent expansion of different types of 
hypergeometric functions with respect to small parameter (in the rest of the paper, we will 
call such expansion as e-expansion) are proposed. They are mainly related to integer values 
of the parameters and/or special values of the argument. The present paper is concerned 
with the Gauss hypergeometric function 



A + ae,B + he 
C + ce 



E 

j=0 



{A + ae)j{B + he)j 
(C + ce), j! 



where {a)j = T{a + j)/T{a) is the Pochhammer symbol, all parameters are real numbers 
and e is a small parameter. Within dimensional regularization [1], the parameter e is related 
with deviation of d-dimensional space-time from its integer value, d = m — 2e. Using the 
well-known representation for the Taylor expansion of the Gamma-function for an integer 
positive number, m > 1 
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or half-integer positive values. 
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where Sk{j) = Y^j=il ^ is the harmonic sum ^, the original hypergeometric function 
with integer and/or half-integer values of parameters {A, B, C} can be written as [2] 
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^The harmonic sums are related with function V-'(^) = 'I^^'^^^i^) ^^'-^ derivatives by means of the 
relation 

r^^'-'Hj) = (-l)'(fc - 1)! [Ck - Skij - 1)] , fc > 1, 

where tp^'^K^) is the fc-th derivative of the T/i-function. In particular, for A: = 1 we have ipiJ) = Si{j — 1) — 7b, 
and 7b is Euler's constant. 
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In this way, the e-expansion of the hypergeometric function p.2|l is reduced to the calculation 
of the multiple sums 
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where mj,nk,luj,Pr are positive integer numbers and \z\ < 1. For z outside of circle of 
convergence it is necessary to transform z in such a way that the original hypergeometric 
functions can be expressed in terms of other hypergeometric functions with convergent series 
[3]. However, very often it is necessary to construct the e-expansion for integrals of the 
following type [4] 
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where all parameters, in general, depend on the e. In this case, it is more convenient to 
construct firstly the ^-expansion of the hypergeometric function and then integrate it with 
the proper kernel. Moreover, by help of the representation 
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the ^-expansion of any generalized hypergeometric function pFp_i can be constructed via 
the ^-expansion of a Gauss one. 

An algorithm for the construction of ^-expansion of hypergeometric functions with integer 
values of parameters has been proposed in [5] and has been generalized recently for rational 
values in [6]. The resulting expansion are expressible in terms of nested sums or multiple 
polylogarithms [7]. 

In contrast to this approach, in paper [8] an alternative algorithm has been invented. 
It is based on construction of e-expansion of a basis of hypergeometric functions ^. For 
hypergeometric functions with integer or half-integer values of parameters, the following 
basis has been analyzed: 
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^The idea of a basis of hypergeometric functions is closely related to idea of master-integrals in high- 
energy perturbative calculations. Using algebraic relations between Feynman integrals derived by help of 
integration by parts approach [9] and/or shifting of the space-time dimension [10], the original set of physical 
amplitudes can be reduced to a restricted set of so-called "master-integrals" . For example, such reduction 
is the necessary step in proving gauge invariance of pole masses (for recent results, see [11]). The algorithm 
of reduction is called the solution of recurrence relations. 
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In this case, the sums (II .3^ are reduced to multiple sums of the following form, 
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where u is, in general, an arbitrary argument (in this particular case it is equal to A'^z) and we 
accept that the notations Sa and Sb will always mean Sa{j — 1) and Sb{2j — 1), respectively. 
For particular values of k, the sums fll.lOj) are called 



harmonic 

1 inverse binomial 
— 1 binomial 



> sums 



The analytical results for harmonic, binomial and inverse binomial sums of different weights 
and depths have been presented in [8,12,13], [8,15], [8,14], respectively. The results of 
the e-expansion are expressed in terms of polylogarithms [16], Nielsen polylogarithms [17] 
or harmonic polylogarithms [18]. The missing part of the approach described in [8] is an 
algorithm for the reduction of original functions to our basis ()1.9j) . However, for all physically 
important cases [19,20], the solutions have been presented. They were derived as the solution 
of recurrence relations for the proper Feynman diagrams [21,22]. 

In this paper we construct an algorithm for the reduction of a Gauss hypergeometric 
function with arbitrary parameters to two Gauss hypergeometric functions with defined 
parameters (reduction to the master-integrals). For integer and half- integer values of pa- 
rameters ^, the e-expansion is constructed up to functions of weight 4. As an illustration 
of the elaborated algorithm, some multiloop scalar integrals are calculated. We like to note 
that for some physically important cases, the proper 5-expansions of Gauss hypergeometric 
functions have been presented in [24-28] . The all-order e-expansion of basis Gauss hyper- 
geometric functions with integer values of parameters is constructed in [29]. The results 
are expressible in terms of multiple polylogarithms. The algebraic Gauss hypergeometric 
functions have been studied in [30]. The numerical approaches are discussed in [31]. 



2 Reduction of Gauss hypergeometric functions 

As is known, for any three contiguous Gauss hypergeometric functions there is a contiguous 
relation, which is a linear relation with coefficients being rational functions in the parameters 
A,B,C and argument z. Using the well know relations [32] 

^Unfortunately, in physical applications other values of parameters also exist [23]. 
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any Gauss hypergeometric function with arbitrary parameters is reduced to the combination 
of eight ones 

■ {a,a + + 1} 
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where a, b, c are some fixed values of parameters. Applying the relations 
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we are able to reduce an original Gauss hypergeometric function to the linear combination 
of two (our basis) 
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These basis functions fj2.7|) are related by a differential identity: 
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In the case when some of the parameters are positive integers (let us put B = m), after 
applying the relation ()2.ip we get one function with the value of one of the parameter equal 
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to unity and some polynomial with respect to z (parameter 5 = 0). In this case, instead 
of relations Eqs. (j2.3j) - ()2.(i|l the following two relations (see [32]) should be used for further 
reduction: 
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In this way, if one of the upper parameters is an integer, then the result of reduction is 
expressible in terms of one Gauss hypergeometric function and a polynomials (the function 
iFq). For case c = b, the relations (j2.3|) and (j2.6|) are useless (0=0). In this case, we should 
apply the Kummer relation ()3.12|) or ()3.15p and reduce this case to the previous one ()2.9|) 
and (jTTn|l : 
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Another algorithm of reduction is described in [33]. 



3 Relations between basis hypergeometric functions 
with integer or half-integer values of parameters 

Let us consider a Gauss hypergeometric functions with integer or half-integer values of e- 
independent parameters. In this case, the set of basis functions consist of the 12 (sixth time 
two) functions. We will call these basis functions as functions of type A, B, C, D, E, F. 
For each type the values of a, 6, c, parameters of our basis ()2.7j) . are presented in Table I: 
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The number of independent basis hypergeometric functions, enumerated in Table I, can be 
reduced by help of the Kummer transformations [34] of variable z. With respect to this 
transformations the functions of type A, B, C, D are transformed into each other. This 
allows us to reduce the number of independent hypergeometric functions. The functions of 
type E, F transform into functions of the same type. Let us illustrate how functions of type 
B, C, D can be expressed in terms of functions of type A. Starting from the relation 
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we express the functions of D-type in terms of functions of A-type 
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Functions of C-type can be written as a linear combination of functions of A-type. Using 
relation 
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As a result, we get the following statement: 

Any functions of type A, B, C, D can be expressed in an algebraic way in terms of just one 
of these types. 

By help of the representation [32], 
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it is possible to find the integral relations between the coefficients of e-expansion of basis 
functions. In our case (integer and half-integer values of parameters) there are integral 
relations between functions of type A and E; B and C; D and F. 

Putting m = 1/2 in the r.h.s. of Eg. ()3.2H1 and using the quadratic transformation 
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we get the representation 
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This representation allows us to find integral relation between coefficients of the e-expansion 
of functions of B and E types; C and E types; and F and A types. 
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4 Laurent expansion of basis functions with integer or 
half-integer values of parameters 

In accordance with the algorithm described in [8], the Laurent expansion of our basis func- 
tions (see Table I in section IH} with respect to parameter e is reduced to the study of multiple 
series of type Using the general expressions given in [8] (see Eqs.(2.30), (2.31), (D.l), 

(D.2) in [8]) we derive that hypergeometric functions of type B and D are expressible in 
terms of multiple binomial sums and hypergeometric functions of type C and E are express- 
ible in terms of multiple harmonic sums. The hypergeometric functions of type A and E 
are expressible in terms of multiple inverse binomial sums. The functions of type F include 
multiple double binomial sums (see Eq. fl4.38j) ). 

Here we present the e-expansion of our basis functions, enumerated in Table I. We restrict 
ourself by constructing the e-expansion up to functions of weight 4. We like to note that 
our expansion is organised in such manner that a term 0{e^) in brackets means functions 
of weight 5. At this order the Nielsen polylogarithms are not enough for constructing the 
^-expansion and new function (one at least) should be introduced. 



4.1 type A 

The e-expansion of functions of type A have been studied in our previous paper [8]. Here, 
we collect the proper results (see Eqs. (2.3), (2.30) and Table I in Appendix C): 
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where 



1 



^2/ [2Li3 {y) - {y)] + 2(ai+a2-/) pLig (-y) - InyU^ {-y)] 

(01 + 02) \TC'y+ (ai + a2-3/)C2ln?/ + (3ai + 3a2-7/)C3| 

|4/(ai + a2-2/)ln|/Si,2(y) 

-2/(ai + a2-/)lni/Si,2(z/^) - 4(ai + a2-/)(ai + a2-2/) lnySi,2(-l/) 
+2/(01 + 02) InyLis (y) + 2(ai + a2)(ai + a2-/) InyLia (-y) 
+2[(ai + a2-/)Li2(-y) + /Li2 {y)f 

-/(ai + a2) In^i/Li2 (y) - (ai + a2)(ai + a2-/) \y? ylji2 {-y) 

+2(ai + a2-3/)C2 /Li2 {y) + (ai + a2-/)Li2 (-y) + -(01 + 02) \t? y 

-^(a? + a2 + «i«2) In^y + Cs Iny 2(ai + a2)^-7/(ai + a2-/) 

+^C4 [l0aia2-30/(ai + a2) + 5(a? + a^)+45/2] I + 0(6^) J , (4.25) 
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(1 + y)^ 1-1/ ci 
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4.2 type B 

For getting the e-expansion of a function of B-type, we will use the representation ()3.19|) 
and ()3.20|) . where r.h.s.'s of the proper equations are given by relations ()4.24|) and ()4.25p . 
In this case, the result of the e-expansion can be written in compact form in terms of the 
variable x related with variable y (j4.26|) as following: 



x = y 



4x 



z^l-l/z ■ 



r, 1-X d 1 + x d 



l + yr^' (1 + x)'' ^' ~ 1+X ' ~dz l-x'^dx 

The form of the variable x automatically follows from our algebraic relations and variable y 
defined by Eq. ()4.2(ij) . Up to order 0{e^) (functions of weight 3) the result of e expansion of 
functions of type B can be cross-checked via multiple binomial sums studied in [8,15]. The 
results of order 0{e^) are new. 



4.3 type C 

In this case we will use the relations ()3.16|) and ()3.17|) for construction of e-expansion. The 
variable y is transformed into variable yc via 

I _ 1- 

yc - yl.^-./d-.) - 1 + ^ • 



(4.28) 
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The form of a new variable yc automatically follows from our algebraic relations and the 
definition of y, Eq. ()4.2(jj) . See also discussion in Appendix D of [8]. For this case, the 
previous results [8] allow us to get expansion only up to functions of weight 2. The next two 
orders of e-expansion are new. 

4.4 type D 

For the e-expansion of D-type functions the relations ()3.13j) and ()3.14|1 are used. In this case 
we have the original conformal variable ()4.2(j|l . Only logarithmic corrections (functions of 
weight 1) are available from the results of [8]. The next three orders of ^-expansion are new. 

4.5 type E 

Functions of type E are algebraically independent from previous ones. Their e-expansion 
can be derived by help of relations presented in [13]. It has the following form (see also 
Eqs.(D.18) in [8]): 



1 + ais, 1 + 



2 + ce 

[ai + a2)c - - 2aia2 Si,2(-2) + 

-c^Lia (2) - i(c - fli - 02)^ ln^(l - z) 
o 



^ ( - ln(l -z)- J 1^2(1 _ ^) + cLi2 (z) 



ai + a2)c — c — aia2 



ln(l - z)U2 (z) 



-sHc 



[ai + 02)0 - - 2aia2 S2,2(-2) + c 



+ (c - ai)(c - a2)(c - ai - 02) 



^ai + a2)c — c — aia2 
1 



ln(l - z)Us (z) 



ln(l- 2)81,2(2) + I ln2(l-2)Li2 (2) 



1 



+ ^(c - ai - a2f \n\l - 2) + c(c - ai - 02)^81,3(2) + e^Lu (2) \ + 0{e*) J ,(4.29) 



l + ce 



1 + aia2e^ Li2 (2) - e< (c - ai - 02)81,2 (-2) + cLis (2) 



+e < c Li4 (2) + (c - fli - 02) 81,3(2) + - [c(c - ai - 02) + 0-102] [Li2 (2)]' 



- [c(c - ai - 02) + 2aia2] 82,2 (-2) \ + C(r 



(4.30) 



In particular, there are two interesting cases: 



as, he 
l + he 



00 i—l 



i-E^*E«'(-^)'"'s.-m( 



(4.31) 



=2 fc=i 
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—e 












— £ 



(4.32) 



1=1 



4.6 type F 

Let us start from the one-fold integral representation 



2F1 
2F1 



l+ce 



l + hiE, 1 + 1)26 

2+ce 



2r(l+ce) 



r(i+6i£)r(i+(c-6i)£ 

2r(2 + c£) 

r(|-i-6i£)r(i+(c-6i)£ 







V2 (sin ^)2''i^(cos 0)2(^-^1)^ 
(l-zsin2^)V2+;.2s ' 

(sin^)^2fei£(pQg^)2{c^i)£ 

(l-^sin2^)3/2+fe2. 



(4.33) 

The finite part of the first integral is equal to complete elliptic integral of the first kind K{k) 
defined as 

f-7r/2 

K{k) ' 



1 — k"^ sin 



2 



1 1 

2' 2 
1 



(4.34) 



The finite part of the second hypergeometric function is expressible in terms of the complete 
elliptic integrals of the first and second kind [36]. Using the relation 

;i+26i£)(l+262e).,,, ^ f l^+he,l+b2e 



-zil-z) 2F1 



2{l+ce) 
= [l + 2(c-6i)e] 2F1 



\ 2+ce 
-^ + bie, ^ + b2e 



2 ' "i^' 2 
l + ce 



+ [{l + 2b2e)z-{l + 2{c-b,)e)] 2F, 



l + bie, \+b2e 
l + ce 



and definition of the complete elliptic integral of the second kind E{k) 

f-7r/2 



E{k) 



we get 



3 3 
2' 2 

2 



1 — k"^ sin^ 



— iF-i 
2 



1 1 

2' 2 
1 



(4.35) 



(4.36) 



(4.37) 



J nz{l — z) 

The next coefficients in the e-expansion of the functions ()4.33|) are related to some gener- 
alization of elliptic functions [37]. In terms of the multiple sums, these new functions are 
related to multiple double sums defined as 



E 



Sai ■ ■ ■ SapSf)^ ■ ■ ■ Sf) 



(4.38) 



In the rest of present paper, we omit functions of type F from our consideration. For the 
definition of coefficients of the e-expansion of functions of type F, the one-fold integral 
representation the ()4.33p or representation of type ()3.21|) and/or ()3.23p should be used. 
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5 Application to Feynman diagrams 



There are several important master-integrals expressible in terms of 2F1 hypergeometric 
functions. This set of integrals includes one-loop propagator type diagram with arbitrary 
values of mass and momentum [26,39]; two-loop bubble integral with an arbitrary values 
of masses [26,39], and one- loop massless vertex diagram with three non-zero external mo- 
menta [26]. For these diagrams, all order e-expansions can be written in terms of Nielsen 
polylogarithms only [26]. Our technology for the expansion of hypergeometric functions 
has been applied also in [8, 15, 38] to more complicated case of generalized hypergeometric 
function, like 3F2 and 4F3. The two-loop propagator type diagram VlOOl and three-loop 
bubble-type diagram E3 (see [26] for details) are expressible in terms of 2-Fi-functions of 
argument 1/4. In this case, the result of expansion can be written in terms of generalized 
log-sine functions and calculated with high accuracy by the help of program LSJK [40]. 

Below we present some diagrams where the results are expressible in terms of Gauss 
hypergeometric functions for an arbitrary set of indices. The proper diagrams are shown 
in Fig. n Diagrams of this type suffer, in general, from irreducible numerator, so that the 



m m 




Figure 1: Diagrams considered in the paper. Bold and thin lines correspond to massive and 
massless propagators, respectively. 



solution of recurrence relations is nontrivial problem (besides one-loop propagator and two- 
loop bubble cases). The solution of recurrence relations for two-loop sunset-type diagram 
was presented in [22,25]. Using the algorithm [10,41], any tensor integral can be presented 
in terms of scalar integrals with the shifted space-time dimension and arbitrary (positive) 
powers of propagators. In scalar integrals of given type, massless subloops can be integrated, 
and the original integrals effectively reduce to more simple integrals with some powers of 
propagators shifted by terms proportional to e. However, for the gauge invariance reason, 
it is desirable to reduce all diagrams to a set of master-integrals before construction of 
£-expansion. 

Using the reduction algorithm described in Sec. 2 it is possible to express the arbitrary 
tensor integral of given type in terms of our basis functions (master-integrals) and integrals of 
more simple structure. The e-expansion of master-integrals can be done by help of relations 
presented in Sec. 4. We would like to note that we are working in Euclidean space-time 
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(p^ — — M^) and use the normalization that each loop is divided by a factor r(l + e). 
5.1 g— loop propagator with q massless line 

Let us consider the g-loop sunset-type propagator with one massive line and q massless ones. 
The result is expressible in terms of Gauss hypergeometric function 

JvAptx, a;2, • • • , a?, /3, , P^) = 

1 r d^hdJ'k2 ■ ■ ■ dJ'kq 



Yi (3-f) J [ki + m?Y[{ki - k2 kg -p)2]"i[A;. 



X2-F1 



ai 



r(c+/3-|g)r(a-f(g-i)) 



a - f (g - l),a+(3-^q 



r(/?)r(f) (3-f) 



(5.39) 



where 



a = ^ a,. 

r=l 



For given type of diagram there are only two nontrivial master-integrals. In the parametriza- 
tion n — 2m — 2e with integer m, the basis is 



l+£{q - 1), l+eq 
2-e 



P F 

2 ' 



e{q - l),eq 
1-e 



p 



(5.40) 



The two-loop case (g = 2) has been considered in [25]. 
5.2 g— loop bubble with g — 1 massless lines 

Let us consider now the g-loop bubble with two massive lines and q — 1 massless lines. The 
result is 



Vi2(q;i, q;2, • • • , Qip, /^i, rn^, M^) = 

ri (3-1) / [A;2 + m2]^ip2_i + M2]/^2[A;2]ai . . . [A;2_2]"-4(A;p - A;i - A;2 VO']"-^ 



d"A;id"A;2 • • • d"A;<, 



^qn/2^^2-^qn/2-Pi-P2-a 



r(/3i)r(/?2)r(f) F'/ (3-f) 

x|f F (a+/32-|(g - 1)) r (^a+(3^+(32-^q) 

' /5i+/32+a-ig,a+/32-i(g-l) 



X2F1 



2^1'^ ' t-'^ 2 ' 
l+/32-f 



15 



+r fe-|)r(a+A-^(g-l))r(a 



X2F1 



a-f(g-2),a+/5i-f(g-l) 



I-P2 + 




n/2-l32 



(5.41) 



where 



a 



q-l 

r=l 



The results of the reduction are expressible in terms of four Gauss hypergeometric functions. 
In the parametrization n = 2m — 2e, where is m is an integer number we get four basis 
functions: 



2F, 



l+e{q-l), 1+eq 
2 + e 

l+e{q-2),l+e{q-l) 
2-e 



iFr 



e{q-l),eq 



l + e 



— 2 ' 2-ri 



^(g-2),£(g- 



-t) ■ (5.42) 



Only for q = 2 (two-loop case) these four hypergeometric functions are expressible in terms 
of one Gauss hypergeometric function and the function iFq, so that only one nontrivial 
master- integral exists. It was calculated in [42]. For q > 2 (3-loop or more) there are four 
independent Gauss hypergeometric functions. As a consequence, there are four nontrivial 
master-integrals for diagrams of this type at 3-loop or more. 



6 Conclusion 

In this paper we have presented the reduction algorithm for Gauss hypergeometric functions 
with arbitrary values of parameters to the two functions ()2.7p with fixed values of parameters, 
which differ from original ones by integers. 

It was shown that the Gauss hypergeometric functions with integer/half- integer values of 
parameters can be divided into 6 types (see Table I). Only three type of them are algebraically 
independent. We have presented the explicit relations which allow us to express the functions 
of type B, C, D in terms of functions of type A (see Eqs. (Km . (imj) . (nTTHll . (ITTTI) . STH^ . 
()3.20|) ). For functions of type A, B, C, D, E the higher-order e-expansion up to functions of 
weight 4 are constructed (see Eqs. ()4.24|1 . ()4.25j) . ()4.29|l . ()4.3()j) ). The result of the expansion 
is expressible in terms of Nielsen polylogarithms only. The e-expansion of function of type 
F is expressible in terms of new functions related to generalizations of elliptic functions. 

As an illustration of the application of our algorithm of reduction, we have considered the 
reduction of g-loop off-shell propagator diagrams with one massive line and q massless lines 
and g-loop bubble with two-massive lines and q — 1 massless lines. We demonstrated, that 
the number of master-integrals for sunset-type diagram beyond one-loop does not depend 
on the number of internal massless lines and it is equal to two. For bubble type-diagram 
beyond two-loop, the number of master-integrals is equal to four. 
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